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ABSTRAGT 

We present a new observational method to evaluate the exponent of the star formation law as initially 
formulated by Schmidt, i.e. the power-law expression assumed to relate the rate of star formation in 
a volume of space to the local total gas volume density present there. 

Total volume densities in the gas clouds surrounding an OB association are determined with a simple 
model which considers the atomic hydrogen as a photodissociation product on the cloud surfaces. The 
photodissociating photon flux incident on the cloud is computed from the far-UV luminosity of the 
OB association and the geometry. As an example, we have applied this "PDR Method" to a sample 
of star-forming regions in M33 using VLA 21-cm data for the Hi and GALEX imagery in the far-UV. 
With these two observables, our approach provides an estimate of the total volume density of hydrogen 
(atomic -I- molecular) in the gas clouds surrounding the young star cluster. 

A graph in logarithmic coordinates of the cluster UV luminosity versus the total density in the 
surrounding gas provides a direct measure of the exponent of the star formation law. However, we 
show that this plot is severely affected by observational selection, which renders large areas of the 
diagram inaccessible to the data. An ordinary least-squares regression fit to a straight line therefore 
gives a strongly biased result. In the present case, the slope of such a fit primarily reflects the 
boundary defined when the 21-cm line becomes optically thick and is no longer a reliable measure 
of the Hi column density. We use a maximum-likelihood statistical approach which can deal with 
truncated and skewed data, and also takes account of the large uncertainties in the total gas densities 
which we derive. The exponent we obtain for the Schmidt law in M33 is 1.4 ± 0.2. 
Subject headings: galaxies: individual (MBS) — ultraviolet: galaxies — galaxies: ISM — ISM: clouds 
— ISM: molecules — ISM: atoms 
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1. INTRODUCTION 

Some 50 years ago, Maarten Schmidt published a pair 
of seminal papers titled "The Rate of Star Formation" 
that discussed the relationship between the interstellar 
gas and the present and past rates of star formation in 
the Galaxy (Schmidt 1959, 1963). This was a topic of 
considerable interest at the time, since radio telescopes 
had begun to map the distribution of atomic hydrogen 
over the Galaxy in some detail. Schmidt postulated a 
simple power-law relationship between the volume den- 
sity of the interstellar gas in a region of the Galaxy and 
the number of stars per unit volume and time formed 
there. He wrote this postulate as: 



SFR( 



(1) 



where SFR is the local star formation rate e.g. in so- 
lar masses per cubic parsec per year, and n is the total 
volume density of the interstellar gas. For the former, 
Schmidt used the Z-distribution and counts of young 
stars. For the latter, he used estimates of the density 
and Z-distribution of atomic hydrogen, while at the same 
time acknowledging that these values might be wrong ow- 
ing to the unknown amount and Z-distribution of molec- 
ular hydrogen which may also be present. With these 
simplifying assumptions, Schmidt found that the expo- 
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nent a in equation [T] was approximately 2 in the local 
Galactic neighborhood. 

The difficulty of obtaining estimates of total gas vol- 
ume density in the interstellar medium (ISM) has pre- 
vented further attempts to obtain estimates of the star 
formation law as initially formulated by Schmidt. In- 
stead, the suggestion was made to relate a star formation 
rate per surface area to the surface density of atomic 
hydrogen as determin ed from the surfa ce brightness in 
the 21-cm Hi line fe.g. lBuat et al.lll989D . While this ap- 
proach necessarily ignores differences in the line-of-sight 
distributions of the gas and young stars, it is straigh- 
forward to implement. For the last ~ 20 years this ap- 
proach has been the mainstay of observational studies 
of the star formation law in disk galaxies at low inclina- 
tions, and the formulation is now called the "Kennicutt- 
Schmidt Law" in recognition of the extensive work in 
this are a by R. Kennicu tt. These studies are summa- 
rized in iKennicuttl (|1998[ ). where it is concluded that the 
exponent a equals 1.4 in an expression for surface den- 
sities of the form T^sfr oc ^gas- I* is then generally 
assumed that the a determined from surface densities is 
the same as the a for volume densities in our equation 
[TJ In more recent studies the 21-cm Hi data has been 
augmented with estimates of the molecular content using 
the surface brightness of the CO(l-O) line as a tracer, but 
the general approach ha s otherwise remained the same. 
iDobbs fc Pringlj ()2009[ ) have recently summarized the 
subject and shown that a simple physical model of star 
formation in the ISM is capable of reproducing a power- 
law form for the Kennicutt-Schmidt law. 
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The existence of a power-law form for the star for- 
mation has been confirmed in many galaxies , althouRh 
a wid e range of values has been reported (|Kennicuttl 
Il997f l. The recent addition of CO as a tracer for the 
molecular gas has not redu ced the scat t er in the re- 
sults. For instance in M33, IHever et all (j2004[) find a 
slope of 1.38 for the molecular gas (based on CO mea- 
surements) but 3.3 for the total gas relation (with a 
monotonically-increasing atomic gas fraction at incr eas- 
ing galactocentric distances) and iVerlev et al.l (|2010D re- 
port values varying from 1.1 to 2.9. A number of concerns 
can be expressed about the meaning of these results given 
the input data. For instance, what are the consequences 
of averaging along the line of sight when we know from 
Galactic studies that the spatial distributions of the trac- 
ers are not the same on scales below a few hundred par- 
sec? And how reliable are the tracers we are using for the 
gas (the beam-smoothed 21-cm and CO(l-O) line surface 
brightnesses) in providing column densities independent 
of local physical conditions and radiative transfer effects? 

In view of such concerns and of the wide range of val- 
ues reported for the exponent in the Kennicutt-Schmidt 
Law, it seems useful to develop other methods for ex- 
amining the quantitative relationship between gas and 
young stars in galaxies. In this paper, we present a new 
method based on treating the Hi found near OB asso- 
ciations as arising in photo-dissociation regions (PDRs) 
which develop on the surfaces of the parent giant molec- 
ular clouds (GMCs) under the action of the far-UV ra- 
diation produced by nearby young stars. This approach, 
which we call the "PDR method", provides us with a 
way of estimating the total gas volume density (atomic 
-I- molecular) in the parent GMCs. The far-UV luminos- 
ity of the nearby OB association is a measure of the star 
formation rate at that location in the galaxy, so that a 
diagram of these luminosities vs. the estimates of volume 
densities in the surrounding clouds (of which there may 
be several) is amenable to interpretation as a star forma- 
tion law in the sense originally described by Schmidt. 

We describe the PDR method in more detail in [J2l 
The data used, and our approaches to dealing with the 
effects of observational selection, are detailed in ij3l We 
estimate the exponent of the star formation law in M33 
in two ways in 21 Our conclusions are summarized in 

m 

2. THE PDR METHOD 

Molecular hydrogen is difficult to observe directly in 
the ISM, owing primarily to its lack of a dipole moment. 
The most common indirect means of inferring its pres- 
ence in the ISM is to observe the low-level rotational 
lines of the carbon monoxide (CO) molecule. The alter- 
native method we present here is based on the physics 
of photodi ssociation o f molec ular hydrogen and the dis- 
covery by lAllen et all ()1986[ ) that this process was re- 
sponsible for the production of atomic hydrogen in the 
promine nt spiral arms of M 83. The method was first fea - 
tured inlA llcn et alJ (|1997D then in [S mith et al.' (2000), 
iHeiner et a l. (20083), and iHeiner et al. (2008b). The 
method provides estimates of the total hydrogen volume 
densities (atomic -|- molecular) in gas clouds located in 
close proximity to regions of recent star formation. These 
clouds may be considered representative of the parent 
GMCs out of which the young stellar association formed. 



In IHeiner et al.l ()2009l ) we presented the first results of 
applying the PDR method to a study of GMCs in M33 
at a linear resolution of 20 pc, the highest available in our 
observations. Here we use the data from a larger survey 
of clouds (Heiner et al., in preparation) at a linear resolu- 
tion of « 80 pc to specifically address the star formation 
law in this galaxy. We begin by describing the geometry 
and reviewing the relevant physics of our model. 

2.1. A simple PDR model 

Clusters of young OB stars are copious producers of 
far-UV (FUV) photons in galaxies, and the FUV lumi- 
nosity of the stellar associat fon is a measure of the local 
star formation rate (see e.g. I Kennicut l 1 1 9981) . The FUV 
photons also radiate into the surrounding interstellar 
medium and dissociate some of the molecular gas found 
in the remaining GMCs. Owing to the high porosity of 
the ISM, the dissociating photons penetrate hundreds of 
parsecs into the int erstellar medium (as confirmed by e.g. 
IHeiner et al.ll2009t ). creating "skins" of atomic hydrogen 
on the surfaces of the GMCs which they encounter; see 
Figure [U As we shall show, the Hi column densities so 
produced on the surfaces of the GMCs can be directly re- 
lated to the total hydrogen volume density in the cloud; 
that gas will be mostly atomic on the cloud surfaces, but 
predominantly molecular deep inside the GMCs. The re- 
sulting total gas volume densities obtained for samples of 
GMCs surrounding a selection of OB associations is then 
combined together with the local FUV-determined star 
formation rate for each association to obtain an estimate 
for the local star formation law in the galaxy. 

A key step in the modeling is to use the physics of 
H2 dissociation by FUV photons in order to relate the 
observed Hi column densities on the cloud surfaces to 
the total gas densities in the clo ud. We use a si mple 
"slab" model des cribed initia lly bv lSternberd (I1988D and 



in more detail by Allen' (2004), with improved coefficients 
provided by Allen et al. (2004). 

The Hi column density in a PDR is calculated with 
the same physics used to determine the excitation o f the 
H2 near-infrared fluorescence lines (lSternberg|ll988D. W e 
use the formulation in Appendix A of lAllen et al. ni200l . 
Briefly, the model is a simple semi-infinite slab geometry 
in statistical equilibrium with FUV radiation incident on 
one side, and a Hi o H2 dissociation-reformation equi- 
librium in the slab on the right side. The solution ap- 
propriate for our present purposes gives the steady state 
Hi column density along a line of sight perpendicular to 
the face of the slab as a function of Go , the inc ident UV 
intensity scaling factor (see Appendix B of A llen et al.l 
120041 for a definition of Go), and the total volume den- 
sity n of H nuclei in the slab. The result is: 
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(2) 



where D is the unattenuated H2 photodissociation rate 
in the average interstellar radiation field, R is the H2 for- 
mation rate coefficient on grain surfaces, a is the effective 
grain absorption cross section per H nucleus in the FUV 
continuum, Gq is the incident UV intensity scaling fac- 
tor, iV(Hi) is the Hi column density, and n is the volume 
density of H nuclei in the slab. This equation has been 
developed using a simplified three-level model ((Sternberg 
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Fig. 1. — A schematic view of a typical region with candidate PDRs. The OB star cluster floods the surrounding highly-porous ISM with 
FUV radiation, and layers of atomic hydrogen are formed on the surfaces of the remaining molecular clouds. 



119881) for the excitation of the H2 molecule. It is appli- 
cable for low-density (n < 10'' cm^^)^ ^^^^ < 590 k), 
isothermal, and static conditions, and neglects contribu- 
tions to A''(Hl) from ion chemistry and direct dissocia- 
tion by cosmic rays. The quantity S' is a dimensionless 
function of the effective grain absorption cross section a, 
the absorption self-shielding function /, and the column 
density of molecular hydrogen A'2 : 



5" = 



N2 



The function S becom es constant for l arge values of N2 
due to self-shielding ()Sternbersj |1988() . Us ing the pa- 
ramet er values in this equation adopted by lAllen et alj 
(|2004D . and writ ing explicit ly the dependence on the 
dust-to-gas ratio 



; mg explicit ly 
(|Allenll2004D . 



we have: 



,,,,, , 7.8 X lO^o , 
^(Hi) = — In 
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(3) 



where A^(Hi) is A^i, the (background subtracted) atomic 
hydrogen column density (in cm^^), S/Sq is the dust-to- 
gas ratio in the ISM expressed in terms of the value Sq in 
the neighborhood of the sun, Gg is the (background sub- 
tracted) incident FUV flux calculated at the affected Hi 
patch, and n = ni+ 2n2 is the total baryonic gas volume 
density (atomic + molecular). This gas will be mostly 
atomic on the surface of the GMC and mostly molecular 
deep inside the cloud. If we can measure A^(Hi), 6/60, 
and Go, equation [3] can be inverted to solve for n; this is 
the essence of our approach to obtaining volume densities 
in interstellar clouds. Although the exponential form of 
the result provides only a noisy estimate of the volume 
density, many measurements of n can be made on GMCs 



in the immediate neighborhood of an OB association, 
thereby improving the overall precision. 

Note that the method we have described does not make 
use of the observed Hi column density to calculate the gas 
volume density from some estimated (but not measured) 
thickness. Such an estimate of Hi volume density ni may, 
however, provide a useful "lower limit" check on the value 
of n = ni -|- 2n2 calculated by the PDR method, since 
the transition region from Hi to H2 is known to be very 
abrupt in most PDRs. 

To summarize, the steps involved in applying the PDR 
method are: 

• Identify the sample of OB star clusters in the 
galaxy around which candidate PDRs can be 
sought. The clusters are knots of bright UV emis- 
sion, and the candidate PDRs are patches of Hi 
surrounding them. 

• Determine the FUV flux of the OB clusters, e.g. on 
GALEX FUV images. The effective wavelength of 
the GALEX FUV band is w 150 nm, and since the 
spectrum is expected to be quite flat in this wave- 
length range the observed GALEX FUV is a rea- 
sonable proxy for dis sociating radiation at « 100 
nm (see lDrain(j| 19780 . 

• Identify the associated Hi patches on the surfaces 
of the surrounding GMCs, measure their surface 
brightnesses on VLA 21-cm images, and convert 
those brightnesses to Hi column densities. Here 
we need to be cognizant of the fact that, with our 
relatively high spatial resolution, we may finally be 
resolving the Hi features, thus becoming sensitive 
to optical depth effects. Indeed, evidence for the 
presence of optically-thick Hi features in our M33 
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data will be presented later. 

• Measure, and as far as possible de-project, the sep- 
aration between the central UV source(s) and the 
Hi patches. This requires that the GALEX and 
VLA-Hi images be accurately aligned. 

• Obtain (or estimate) dust-to-gas ratios in the gas, 
ideally local to the candidate PDR. 

This procedure allows us then to calculate the total hy- 
drogen volume densities as local spot measurements at 
the location of candidate PDRs. 

3. DATA 

The Hi observations of M33 used in this paper were 
provided by David Thilker and Rob Braun (2007, pri- 
vate communication) and are presented and discussed in 
Heiner et al. (in preparation) . We used the PDR method 
on these data to calculate the total hydrogen volume den- 
sities used here. The linear resolution of the Hi data is 
« 80 pc. The local star formation rate at the location of 
the candidate PDRs is estimated using the far-UV lumi- 
nosity of the parent OB associations from the GALEX 
data, assuming that the star formation rate is directly 
proportional to the lum inosity at t he GALEX fa r-UV 
wavelength (see lKennicu tt 1998; Ma dau et aLlll998[ )). In 
that case, using the far-UV luminosity results in the same 
power law slope when relating Ljjv to n. In the absence 
of region-specific extinction measurements, we have not 
made any corrections for extinction; a global extinction 
would merely shift all log-luminosities equally, and there- 
fore would not influence the calculated value of the ex- 
ponent in the star formation law. The volume Schmidt 
Law correlation we seek is of the form: 

logL;7y oc alogn. (4) 
4. RESULTS 

The computed values of Luv and n from Heiner et 
al. (in preparation) are shown in Figure [2] in log coordi- 
nates, along with several fits to straight lines. This figure 
shows a rough correlation, as might be expected from a 
star formation law of the form of equation 31 The circles 
correspond to regions inside i?25, whereas the triangles 
are candidate PDRs outside of i?25- We have made no 
attempt to fit these regions separately. The straightfor- 
ward ordinary least-squares (OLS) fit to the data rep- 
resented by the dash-dotted line in the left panel yields 
a = 0.6. 

4.1. Selection effects 

The data plotted in Figure [2] clearly avoid large regions 
of the graph, as follows: 

4.1.1. Low values of Ljjv 

Our survey of star clusters in VI33 cuts off below a 
limit of Luv ~ 2 x 10^^ ergs/sec/ A). VIost of these faint 
clusters (which contain few O stars) are found outside the 
i?25 isophotal radius in regions of low confusion; here the 
lower limit is set by the sensitivity of the GALEX data. 
Such faint clusters also exist inside R25 (along with a 
profusion of widely-scattered FUV-producing B stars), 
but confusion prevents their accurate tally. However, 



complete counts are not required for our study, since we 
subtract away the mean background levels on the FUV 
and Hi images in order to establish the excess Ljjv of 
the cluster and the excess Hi which the cluster's FUV 
photons produce on the surfaces of nearby GMCs. In this 
paper we have therefore not pursued FUV-faint objects 
inside the main disk of the galax}0. 

4.1.2. High values of Luv 

At the bright end of the FUV luminosity we see that 
there are no OB associations more luminous than NGC 
604 at Luv w 5 x 10'^^, the equivalent of several do zens of 
04 stars (as inferred from iBruhweiler et al. I [20031 ). This 
limit is apparently set by the physics of "remaining" steps 
in the star formation process in M33, which we do not 
consider here (nor even pretend to understand). 

4.1.3. Region II 

The sloping dashed line (partially) bounding the empty 
triangular region labelled II in the lower right corner of 
the left panel of Figure [5] indicates an Hi column density 
lower limit of « 6 x 10^^ cm~^, at a typical dust-to-gas 
ratio of around 0.33. The available data progressively 
disappears as we approach this line, which marks the 
approximate sensitivity limit of our VLA-}li data. This 
is a common (but largely benign) observational selection 
effect; the computed values of n are more noisy here, but 
they are not strongly biased. A small bias could creep 
in if the slope of the line reflecting the selection effect 
were significantly different from the slope of the actual 
star formation relation, and at the same time there is 
an abundance of data points near this selection limit. 
Although the former is likely to be true, the latter is 
most assuredly not. 

4.1.4. Region I 

The other empty triangular region labeled I in the up- 
per left corner of the left panel of Figure [5] appears to 
arise from a more interesting (and also a more sinister) 
observational selection effect. The sloping dashed bound- 
ary line here is defined by an Hi column density upper 
limit of 5 X 10^^ cm~^ at a characteristic dust-to-gas ra- 
tio; higher column density estimates are apparently very 
rare. We suggest that this is a consequence of the 21-cm 
line becoming progressively more optically thick at about 
this limiting value, thereby underestimating the true Hi 
column density and hence, following equation |3l over- 
estimating the volume densities n. Points can lie some- 
what to the left of this boundary owing, for example, to 
different dust-to-gas ratios compared to the typical value 
of 0.25 that we adopted. The different typical dust-to-gas 
ratio from the one we used to delineate region II (0.33) is 
a reflection of the fact that the higher column densities 
occur in the inner regions of M33, where the metallicity 
is higher as well. On the other hand, the lower column 
densities are measured mostly in the outer regions, where 
the metallicity is systematically lower. This selection ef- 
fect is an optical depth limit, appearing as a consequence 

A quantitative study of these "backgrounds" is nevertheless 
highly relevant to the interesting question of what fraction of the 
total Hi content of the entire galaxy arises through photodissocia- 
tion, i.e., what fraction of the Hi present is not "primordial", but 
rather a product of the star formation process. 
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Fig. 2. — The far-UV luminosity of our candidate PDRs is plotted against the total hydrogen volume density. Left panel: The selection 
effects are indicated with dashed lines bounding areas where the Hi column density upper limit (1) and the til column density lower limit 
(If) are found. The ordinary least-squares (OLS) fit is represented by the dashed-dotted line, which is Luy oc rfi'^. Right panel: The 
dashed line labeled "Hist" graphs Lijy oc n^'^, obtained using histograms of the distribution of n; see text for further explanation. The 
solid line labeled "MCMC" corresponds to Luy oc n^ *, from the maximum likelihood estimate. This last value of 1.4 (with an estimated 
uncertainty of 0.2) is our best estimate of the exponent in the Schmidt Law for star formation in M33. 



of the improved spatial resolution of the present data set 
that leads to a (partial) resolution of the M33 GMCs in 
Hi. 

4.2. A simple correction 

The appearance of an optical depth limit to the 21-cm 
surface brightness leads to a serious bias in determining 
the true exponent of the star formation law from the 
observations. The dash-dot OLS fit line in the left panel 
of Figure [2] clearly shows a slope (of 0.6) which is too low, 
a consequence of biasing the calculated gas density to 
values that are artificially high. To explore this selection 
effect further, we divided the Luy data into six (log) 
luminosity ranges, and computed the histograms of the 
logn values in these six horizontal strips; the results are 
shown in Figure [H These distributions appear to be 
approximately log-normal, but are also clearly biased by 
the selection effect. As a first attempt to correct for this 
bias, we identified the value of log n corresponding to the 
peak of the histogram in each strip (shaded darker gray 
in the figure) and fitted an OLS line to these six values. 
The exponent found by this (oversimplified) method is 
1.2 (graphed as the dashed line labeled "Hist" in the right 
panel of figure [2]) . This value must still be too small, 
since we have assumed that the histograms of figure [3] 
are merely cut off to varying degrees on their left sides, 
rather than being cut off and skewed towards their right 
sides in some nonlinear way by the bias. Removing the 
outer bins does not change the slope. 

4.3. A regression analysis 



We have shown that selection effects truncate and skew 
the observed distribution of cloud densities, biasing the 
value of the exponent obtained by a simple least-squares 
fit to the star formation law to values that are too low. 
These observational limits cause selection effects similar 
to the well-known "Malmquist bias". In order to ob- 
tain a more statistically-sound value for the exponent a, 
a regression analysis based on the use of "Monte Carlo 
Markov Chains" (MCMC) can be used on the cloud den- 
sity and the UV luminosity data (Equation [3]) , taking 
into account both the bias and the (relatively large) un- 
certainties in our cloud densit y measurements. This form 
of Bayesian inference (see e.g. lTolstov fc Sahal ll996^. uses 
a priori assumptions of the data, that are adjusted by 
the measured data. It provides a statistical formalism to 
obtain the true exponent based on simulations of unbi- 
ased data, where the "simple approach" in our previous 
section provided an intuitive way to estimate the miss- 
ing data with similar assumptions. In both cases the 
assumption is the log-normal distribution of the values 
of n in limited ranges of the UV lu minosit y. 

The fuU procedure is described in iKellvl (|2007l ). where 
it is shown how this approach (referred to as a maxi- 
mum likelihood estimate MLE using a Gaussian mixture 
model) is robust in the presence of censored data and 
large measurement errors. We used the IDL routines de- 
scribed in that paper to perform our analysis. The result 
of the MLE is an a posteriori median estimate, which 
means that the slope (and intersect) resulting from this 
method are direct measurements of their true values. 

The IDL routine constructs a Monte Carlo Markov 
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Fig. 3. — The frequency of occurrence of log(n) in each range of the (log) far-UV luminosity is plotted here. The bin with the highest 
frequency is dark gray. Selection effects cut off the data points on the left side and skew the histograms towards the right. In order to 
illustrate a (log-)normal distribution which might be expected from stochastic noise, dashed lines represent Gaussian distributions of an 
arbitrary height and width centered at the highest frequency bars. 



Chain (| Jewell et al.l[200l review the formalism in their 
Appendices), where new data points are simulated us- 
ing the "Metropolis - Hastings" algorithm. This partic- 
ular algorithm is appropriate when errors dominate, and 
when the variables may not be independent, as is the case 
here. The regression parameters (generally) converge af- 
ter several thousand iterations to provide the parameters 
drawn from the posterior distribution. It is important to 
realize that the value for the exponent in the star for- 
mation law obtained in this way is a random draw from 
the true distribution of all exponents consistent with the 
noisy and censored input data. The method can there- 
fore also be used repeatedly to provide an estimate of the 
stochastic error in the exponent. 

To use the MCMC method as presented by iKellvl 
(|2007l ). we need to adopt reasonable estimates for the 
(la) errors in log(n) and log(L;7y). This is essential 
for the algorithm to estimate the true distribution of n 
and . The crude but intuitive attempt at correcting 
for the selection effects described in the previous section 
shows that the cloud densities (ignoring the truncation 
and skewing) roughly follow a log-normal distribution. 
Also, the errors in n dominate the correlation between 



log(n) and log(L[/y). We take absolute uncertainties in 
the values of log(n) and \o%[Luv) of 0.5 and 0.2 respec- 
tively, conservative values reflecting the relatively large 
error in n and the smaller error in L, where L is directly 
based on the measured UV flujl^- We then run the regres- 
sion routines repeatedly. The resulting slopes are shown 
in Figure [H The distribution of the values of the slope 
is non- normal, but the average value is 1.4 with a range 
of ±0.2. Note that the MCMC method is invariant to 
swapping the coordinate axes. It is also robust against 
small variations in the estimated uncertainties for n and 
and setting these uncertainties to nearly brings 
back the OLS bias. The solid line with slope 1.4 labeled 
MCMC in the right panel of Figure[2]reflects the outcome 
of our Bayesian approach to correct for our selection ef- 
fects, and provides our best estimate of the exponent in 
the Schmidt Law of star formation. 

5. DISCUSSION AND CONCLUSIONS 

The new observational method to estimate the slope 
of the star formation law that we present here can be 

^ Our forthcoming paper on the input data set will have the 
values of observed UV fluxes. 
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Fig. 4. — The distribution of slopes obtained by applying the 
MCMC method 40 times is shown in this histogram. The most 
common value is 1.4, followed by 1.2, while a range of values be- 
tween 1.2 and 1.6 is obtained. 

further improved in various ways. While some selection 
effects may have been dealt with acceptably, others re- 
main, and all such effects need to be considered carefully 
as we have attempted to do here. 

1. The method yields cloud densities with a relatively 
high uncertainty (we used a conservative estimate 
of 50%). This uncertainty could be reduced by im- 
proving the quality of the data; however, the spread 
in cloud densities may very well be intrinsic, similar 
to the scatter in the local metallicities measured at 
different places in the disks of nearby galaxies (e.g. 
iRosolowskv fc Sin^l2008[ ). 

2. The star formation rate is derived from the UV 
fluxes of OB associations, and these may be af- 
fected by extinction to various degrees which we 
have ignored here. A detailed extinction map of 
M33 would be helpful, as only a spatially-varying 
extinction correction would influence our determi- 
nation of the exponent in the star formation law. 

3. A better calibration of the observed UV flux to 
the actual star formation rate could be carried out. 
However, IMadau et al.l (jl998[) found that the star 
formation rate is directly proportional to the FUV 
luminosity at 150 nm, which is very close to the 
effective wavelength of the GALEX far-UV data, 
so it is unlikely that further precision here would 
change our results. 

4. A correction for the effects of high optical depth in 
the 21-cm line could be developed in order to obtain 
more accurate estimates of the Hi column density. 
While such a correction would surely help, it is un- 
fortunately not clear how one might determine it. 



Abramova, O. V., & Zasov, A. V. 2008, Astronomy Reports, 52, 
257 

Allen, R. J. 2004, in Astrophysics and Space Science Library, Vol. 
319, Penetrating Bars Through Masks of Cosmic Dust, ed. 
D. L. Block, I. Puerari, K. C. Freeman, R. Groess, & E. K. 
Block, 731 



5. Lower atomic hydrogen column densities can only 
be used if they can be distinguished clearly from the 
general Hi background. More sensitivity and high 
spatial resolution would be required here, although 
we have argued that the biasing effect of limited Hi 
sensitivity on the power law slope is minor. 
Unique to our approach is the use of cloud volume 
densities, not surface densities, and the enhanced sen- 
sitivity of the PDR method to lower gas densities (less 
than a few hundred baryons cm~^) when compared to 
the higher-density sensitivity of the CO(l-O) estimates. 
T o our knowledge, the only e xception to this is the work 
of lAbramova fc Zasovl ()2008f ): these authors specifically 
aimed to explore a volume density Schmidt Law. How- 
ever, they could make only rough (CO-based) gas vol- 
ume density estimates, based on column density mea- 
surements and an (uncertain) theoretical model of the 
thickness of galactic disks. 

To conclude, we have carried out the first direct de- 
termination of the exponent in the star formation law as 
initially formulated by Schmidt over 50 years ago. The 
star formation rate has been estimated using the far-UV 
luminosities of OB associations in M33. The total gas 
volume density (atomic -|- molecular) has been obtained 
using a method which regards the atomic hydrogen as the 
dissociated "skins" of the molecular clouds out of which 
the young stars have formed. We have used a simple 
one-dimensional slab model for the dissociated layer of 
Hi on the cloud surfaces. We also show how observational 
selection leads to a serious bias in the determination of 
the exponent of the star formation law. We have used a 
"maximum likelihood" technique to account for this se- 
lection bias. Our result for the slope of the star formation 
law in M33 is 1.4 ± 0.2. We note that in this paper we 
have addressed only the question of the value of the ex- 
ponent in the Schmidt Law, and not the normalization. 
Application of our method to a larger sample of galaxies 
will provide additional information on this normalization 
and its universality. 

Finally we point out that the value we have obtained 
for the exponent in the Schmidt Law for star formation 
is close to that used in many recent successful numeri- 
cal simulations of galaxy formation and evolution. This 
correspondence lends further credence to our view that 
the Hi in galaxy disks is a product of the star formation 
process, and that this process can be usefully analysed 
using the simple physical model described in this paper. 
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